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SUM  MAR  Y 

A  further  modification  and  enhancement  of  the  doublet  lattice  computational  proce¬ 
dure  in  use  at  ARL  is  described.  The  theoretical  basis  for  the  new  procedures  used  to 
integrate  the  unsteady  Kernel  function  is  covered  in  detail.  A  number  of  cases  involving 
nonplanar  and  nonparallel  lifting  surface  combinations  are  studied,  and  the  results  are  com¬ 
pared  with  those  of  other  workers.  These  show  generally  good  agreement,  although  the 
predictions  of  the  magnitudes  of  the  forces  acting  on  a  lipstore  and  an  underwing  store 
do  not  compare  as  favourably.  The  present  method  has  been  implemented  at  ARL  into  the 
FORTRAN  77  program  called  DO U I  AT. 
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A ,  B  points  at  either  end  of  the  doublet  line  on  the  sending  panel 

Ak  nondimensional  area  of  the  fc’th  panel  defined  in  Equations  6-25  and  6  26 

C  control  point  on  receiving  panel 
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c  mean  geometric  chord  of  NLR  wing;  c  =  0-4183  m 
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Equation  2-2 

Dil\  D{1)  planar  and  nonplanar  parts  of  £>„  defined  in  Equations  4- 12  and  4-13 

rs  r* 

F  frequency  of  oscillation  in  Hz 

/  nondimensional  mode  shape 

G111,  G(J|  numerators  of  the  planar  and  nonplanar  parts  of  the  Kernel  function  defined  in 
Equations  A  -4  and  A  •  5 

h  vertical  separation  between  wing  and  tail  defined  in  Section  61 

i  s-\ 

i,  j,  k  orthogonal  unit  vectors  in  the  v.  y  and  r  directions 

K  Kernel  function  defined  in  Equation  3- 1 

Afrj  Kernel  function  relating  the  induced  normalwash  at  a  receiving  point  r  to  the 

pressure  at  a  sending  point  s 


Klt  K2  elements  of  the  planar  and  nonplanar  parts  of  the  Kernel  function  defined  in 

Equations  3  -4  and  3  5 

k  nondimensional  frequency  parameter,  wl/U 

I  reference  length 

/,  quarter  chord  line  of  s’th  panel  defined  in  Equation  2-4 

M  free  stream  Mach  number 

M  unsteady  pitching  moment  defined  in  Equation  6  -24 

A/j  unsteady  spanwise  sectional  pitching  moment  defined  in  Equation  6  16 

N  unsteady  yawing  moment  defined  in  Equation  6-33 

n  number  of  straight  line  segments  into  which  the  interval  /„  is  divided 

p  pressure  defined  in  Equation  2- 1 

Q  generalised  force  coefficient 

r  control  point  on  the  receiving  panel  located  at  the  3  4-chord  point  and  illustrated 

in  Figure  4 

rAB(<r)  position  vector  of  points  along  the  doublet  line  of  a  panel  as  a  function  of  non- 

dimensional  position  parameter  a  defined  in  Equation  5  I 

5  semi-span 

7,.  T:  elements  of  the  planar  and  nonplanar  parts  of  the  Kernel  function  defined  in 

Equations  3-2  and  3-3 

T2*  element  of  the  nonplanar  part  of  the  Kernel  function  defined  in  Equation  4-9 

r  time 

U  free  stream  velocity 

ff'7n,  ff'72’  weights  used  in  integration  of  the  planar  and  nonplanar  parts  of  the  Kernel 
function  when  r2  >  0,  defined  in  Equations  A  37  and  A  3X 

m  dimensional  normalwash  defined  in  Equation  2- 1 

.VAB|u).  >AB(a).  components  in  the  i.  j.  k  directions  of  the  position  vector  r.Mt;«)  defined  in 

Zab(<t)  Equation  5  I 

v.  v.  ;  orthogonal  coordinates  of  a  point  defined  using  the  global  (  V.  )',  /)  coordinate 

system  of  Figure  I 

vBM  dimensional  v-ordinate  of  the  balance  centre  defined  in  Equations  6  -26  and  6-35 


•v, 


dimensional  v-ordinate  of  the  lift  point  on  the  A  th  panel  defined  in  Equations 
6  26  and  6  35 


dimensional  x-ordinate  of  the  1 /4-chord  point  of  the  /th  section  defined  in 
Equation  6  18 


unsteady  side  force  defined  in  Equation  6-32 

unsteady  normal  force  defined  in  Equations  6-23  and  6-30 


nondimensional  normalwash  defined  in  Equation  2-3 

parameter  used  as  a  criterion  for  selecting  the  appropriate  integration  scheme, 
defined  in  Equations  5  - 1 5—5  - 1 7 

v'l  -  M1 

sweep  angle  of  quarter  chord  line  of  the  r'th  sending  panel  defined  in  Equation  2-4 

dihedral  angles  of  the  receiving  and  sending  points  defined  in  Equations  3-2 
and  3-3 

angular  rotation  in  torsion 

nondimensional  pressure  defined  in  Equation  2-5 
density  of  air  in  the  free  stream 

parametric  position  coordinate  defined  in  Equation  5- 1 
circular  frequency  of  oscillation 


Subscripts  and  superscripts 


denotes  planar  part  of  variable  or  function 
denotes  nonplanar  part  of  variable  or  function 
denotes  value  at  point  A  on  doublet  line 
denotes  value  at  point  B  on  doublet  line 

denotes  quantity  relating  to  a  panel  in  thc/'th  row  of  the  /‘th  streamvvise  column 
of  panels 

denotes  receiving  panel  or  point 


denotes  sending  panel  or  point 


denotes  a  quantity  relating  to  the  steady  Kernel  function. 


I.  INTRODUCTION 


A  previous  report  by  Waldman  [1]  describes  a  FORTRAN  program  for  calculating  general¬ 
ised  air  forces  on  arbitrary  combinations  of  lifting  surfaces.  The  method  used  therein  is  a 
variation  of  the  doublet  lattice  method  of  Albano  and  Rodden  [2].  This  report  describes 
further  modifications  and  enhancements  to  the  integration  procedures  used. 

The  initial  formulation  of  the  doublet  lattice  method  [2]  is  adequate  for  calculating  inter¬ 
ference  effects  on  general  nonplanar  lifting  surface  configurations;  however,  some  computational 
difficulties  have  occurred  for  nearly  coplanar  wing/horizontal  tail  combinations  [3,  4.  5],  The 
formulation  of  the  doublet  lattice  method  used  in  [1  ]  is  able  to  deal  with  nearly  coplanar  surfaces, 
up  to  the  point  where  the  vertical  separation  is  so  small  as  to  be  negligible.  Although  the  method 
yields  accurate  results,  it  is  computationally  expensive,  especially  for  nearly  coplanar  surfaces, 
since  the  number  of  integration  points  needed  is  large. 

A  significantly  improved  method  for  carrying  out  the  required  integrations  has  been  devel¬ 
oped  and  implemented.  This  uses  some  features  of  the  approach  proposed  by  Giesing  et  at.  [4] 
and  extends  the  techniques  found  in  [)]. 


2.  THE  DOUBLET  LATTICE  METHOD 

The  doublet  lattice  method  is  a  panel  method  for  the  solution  of  the  oscillatory  subsonic 
pressure-normalwash  integral  equation  for  multiple  interfering  surfaces 
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where  n(.v„  :,)  e"'"  is  the  induced  oscillatory  normalwash.  />(.v,.  i;.  :J  c1"'  is  the  pressure 

distribution  over  all  lifting  surfaces.  L.S.,  and  K  is  the  subsonic  nonplanar  Kernel  function. 
The  coordinates  of  the  sending  and  receiving  points  are  given  bv  (v,.  v„.  :J  and  (,vr.  vr.  The 

f 

symbol  «  indicates  integration  in  the  sense  of  Maneler  [6], 


In  the  doublet  lattice  method  the  lifting  surfaces  are  divided  into  small  trape/ial  panels, 
as  shown  in  Figure  1.  By  assuming  that  the  unknown  pressure  p  is  uniform  over  a  panel.  Albano 
and  Rodden  [2]  have  shown  that  the  integral  equation.  Equation  21.  reduces  to  a  set  of  linear 
simultaneous  equations.  These  may  be  written  in  matrix  form  as 


I*,!  =  [/\J  IzJ 
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where 
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4.  BEHAVIOUR  OF  THE  KERNEL  FUNCTION 

The  general  form  of  the  nonplanar  Kernel  function.  A,  was  given  in  Equation  31  as: 

K  =  (A,  7',  f  A\  7, )  /  ’  4-1 

The  formal  basis  of  the  doublet  lattice  method  is  to  approximate  the  numerator  of  A  with  a 
polynomial  and  integrate  Equation  2  3 analytically  (2 ).  The  integration  used  to  obtain  the  normal- 
wash  influence  coefficients.  Dr>.  works  well  for  all  cases,  planar  and  nonplanar.  with  the  exception 
of  the  nearly  coplanar  case  [3.  4.  5J. 

Consider  the  case  where  the  receiving  point  is  downstream  of  the  sending  panel,  where  A 
defines  the  downstream  direction  in  Figure  I.  When  there  exists  a  small  normal  separation 
between  the  receiving  point  and  the  plane  of  the  sending  panel,  then  the  numerator  of  the  Kernel 
function  has  large  variations  along  the  interval  /.  in  the  line  integral 


4  2 


The  typical  behaviour  of  the  numerator  of  the  steady  Kernel  line  non  is  illustrated  m 
Figure  2.  Values  of  the  numerator.  (A  •  A:'"'/.i.  are  plotted  against  the  spanwise 

variable  r,  in  the  vicinity  of  c,  =  0  for  three  values  of  vena  I  distance  and  a  chord  wise 
distance  v,  =  0-5.  Both  the  sending  and  receiving  panels  he  in  planes  parallel  to  the  horizontal 
plane.  Note  that  the  curves  are  asymptotic  to  the  value  (A,"  /,  •  A  /  i  I  n  lor  values 
of  g,  -  ±  x .  although  this  is  not  apparent  from  Figure  2.  (hcsing  ii  <;/  |4J  li.m  found  'fiat 
when  the  length  /,  of  the  line  integral  is  large  compared  with  <  i  '  the  sanations  ot 

(A, 1,1  T,  +  K2"‘  T2)  across  the  panel  are  verv  large  I  hereton  ,i  'u.'ii.l  .'ido  polv  nominal 
fit  is  inadequate. 

t 


Expressions  for  the  various  parts  of  the  unsteady  Kernel  function  valid  for  small  values 


of  yt,  and  downstream  of  the  sending  point  are  [4] 

Kl  -* 2  4-3 

K2  -*  -4  4-4 

f,  =  cos  (y,  -  y5)  4-5 

T2  =  (z,  cos  y,  -  sin  y,)(z,  cos  ys  -  y,  sin  ys)/r2  4-6 


where  r2(  =  y,2  +  z,2)  -» 0  and  x,  >  0.  Using  Equations  4-3  to  4-6  in  the  numerator  gives 


-  m mt  i 

lim  e~~ ~  {X,  T,  +  K2T2} 

r2  -»■  0 


-  itux  i 

e  ~  {2cos(yr-ysH4(z,  cos  yr— y,  sin  yr)(z,  cos  ys— >>,  sin  ys);r2J  4-7 


The  term  that  requires  the  most  care  during  integration  arises  from  the  nonplanar  term  and  is 
the  one  divided  by  r 2  and  involves  T2.  A  plot  of  T2  as  a  function  of  y,  is  given  in  Figure  3  for 
y,  =  0°,  yr  =  45°  and  two  values  of  z,.  It  is  clear  that  a  low-order  polynomial  fit  will  not  give 
accurate  results  when  r/l%  is  small. 

One  solution  to  this  problem  is  to  consider  the  nonplanar  terms  separately  from  the  planar. 
The  planar  terms  vary  as  1/r2  but  the  nonplanar  terms  vary  as  l.r4.  hence  the  Kernel  function 
can  be  written  as  [4] 


-  iiux  t 

K  =  e  22  (K|  T,  r2  +  K2  T2*  r*)  4-8 

where  T 2*  =  r2  T2  =  (z,  cos  y,  —  y,  sin  •/,)(;,  cos  •/,  —  r,  sin  ys)  4-9 

Equation  4-8  may  be  rewritten  as 


K  =  e  f  K,  T,  r  +  e 


4-10 


and  the  numerators  of  each  of  the  terms  may  now  be  approximated  with  a  low-order  poly¬ 
nomial  since  they  are  slowly  varying  functions  over  the  line  interval  /,. 

Hence,  Equation  4-2  may  be  written  as 
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By  a  suitable  choice  of  integration  scheme.it  is  possible  to  analytically  incorporate  the  1/r2 
and  1/r4  behaviour  into  the  integration  procedure.  This  has  the  effect  of  minimising  compu¬ 
tational  errors  involved  in  evaluating  the  integrals  over  regions  where  r  is  small  or  becomes 
equal  to  zero. 


5.  INTEGRATION  OF  THE  KERNEL  FUNCTION 
S.l  Geometric  Considerations 

In  calculating  the  aerodynamic  influence  coefficients  defined  by  Equation  2  -4.  let  us  consider 
the  interaction  of  a  receiving  panel  r  and  a  sending  panel  s.  as  shown  in  Figure  4. 

The  position  vector  corresponding  to  points  along  the  doublet  line  may  be  written  in  terms 
of  the  parametric  variable  a  as  follows 


=  Tab(ct)  i  + 

Eab(">  j  +  AbI^)  k  0  "  1 

5-1 

where 

Tab(<7) 

=  vA  +  (  vB  - 

vA)<r 

5-2 

=  y  \  +  (v’b  — 

>'a>  o 

5-3 

=  -A  +  <-B  — 

-a)  a 

5-4 

From  the  definitions  of  and  z,  in  Equation  3  -6  it  is  clear  that  when  evaluating  the 

line  integral  of  Equation  2-4  the  variables  .v,.y,,  and  z,  may  be  written  in  the  following  form 


.v,  =  a,  +  bxo  v,  =  a,  +  bya 

-i  =  ",  +  V 

5-5 

w  here 

a.  =  v,  —  vA  a.  =  y,  —  t  A 

",  =  ~  -  A 

5-6 

=  va  —  VB  =  3  A  —  .'ll 

=  -A  —  -B 

5-7 

In  the  equation 

for  the  Kernel  function  the  variable  r: 

is  present.  The  equation 

for  r  may 

be  written  as 

r~  a,  o'  +  b,  o  +  c, 

5  -K 

where 

a,  --  by  +  by 

5-9 

h,  =  2[a,  />.  v  a,  b,) 

5  10 

c,  =  ay  -  ay 

511 

If  the  minimum 

value  of  r:  tor  a  given  panel  combination  i->  denoted  bv  r*m, 

in.  and  the 

location  of  this  minimum  bv  then  we  have  that 


5  i : 


4  a, 


5 


5  13 


The  span  of  a  given  sending  panel,  A s,  may  be  calculated  from 

A  s2  =  b,2  +  bx2  5  14 

Due  to  the  organisation  of  panels  into  streamwise  columns,  omjn,  r2mi„  and  As  need  only  be 
calculated  for  combinations  of  columns,  rather  than  individual  combinations  of  panels.  This 
strategy  reduces  the  number  of  times  these  three  parameters  are  calculated. 


5.2  Basis  for  Choice  of  Integration  Procedure 

A  decision  which  determines  the  integration  procedure  to  be  used  is  made  on  the  basis 
of  the  values  rmjn,  omia  and  As  that  have  been  calculated  for  a  given  combination  of  sending  and 
receiving  columns  of  panels.  One  of  the  following  conditions  can  occur  for  such  a  combination 


(•)  rmJAs  >  e 

-x  >  aml„  >  +  x 

5-15 

(ii)  0<  rmJAs<  e 

<rm,„  <0  or  <rmi„  >  I 

516 

(iii)  0  <  rmin/Ar  <  e 

0  <  omi„  <  I 

5-17 

where  e  is  some  small  number  (typically,  6  =  0  01  in  the  program).  The  ratio  rmin/  As  is  a  non- 
dimensional  parameter  that  compares  the  smallest  value  of  r  obtained  for  a  given  combination 
of  sending  and  receiving  panels  with  the  span  of  the  sending  column  of  panels.  For  compu¬ 
tational  purposes,  the  singularity  r  =  0  is  assumed  to  occur  when  rmin  lies  in  the  range  defined 
by  0<  rmJAs<,  e. 

Condition  (i)  corresponds  to  the  case  where  no  singularity  occurs  for  the  combination  of 
sending  and  receiving  panels  (or  columns)  being  considered.  The  integration  procedure  used  is 
described  in  Appendix  A  for  the  case  where  the  interval  /,  is  subdivided  into  n  smaller  sub¬ 
intervals.  In  the  program  the  integration  for  condition  (i)  is  evaluated  with  n  =  1. 

Condition  (ii)  corresponds  to  the  case  where  a  singularity  occurs,  but  it  lies  outside  the 
interval  of  integration  /,.  The  integration  procedure  described  in  Appendix  A  is  applied  w  ith  //  =  1. 

Condition  (iii)  corresponds  to  the  case  where  the  singularity  occurs  within  the  interval  of 
integration  /,.  The  singularity  is  assumed  to  occur  in  the  centre  of  the  interval.  This  requirement 
can  be  met  by  aligning  all  columns  of  panels  in  streamwise  strips.  The  integration  procedure 
used  is  described  in  Appendix  B  for  the  case  where  the  interval  /„  is  subdivided  into  (2n  —  1) 
straight  line  segments.  In  the  program  n  =  4  and  the  singularity  falls  in  the  centre  of  the  4th 
subinterval.  The  width  of  the  fourth  subinterval  is  !s  n.  The  remaining  six  subinlervals  are  placed 
symmetrically  about  the  segment  containing  the  singularity,  three  on  either  side,  and  the  length 
of  each  of  these  subintervals  is  IJ(2n).  The  choice  of  n  =  4  follows  the  technique  used  by  Farrell 
[5],  By  changing  a  PARAMETER  statement  in  the  FORTRAN  source  code  (IR2EQ0  =  2 n. 
n  =  1,  2,  .  .  .)  it  is  possible  to  vary  the  integration  scheme  used  to  deal  with  the  singularity. 

In  order  to  ensure  that  a  high  degree  of  numerical  accuracy  is  maintained  in  all  cases,  it 
is  necessary  to  choose  an  appropriate  value  of  the  parameter  < .  used  previously  in  Equations 
5-15  to  5-17.  If  the  chosen  value  of  e  is  too  large  then  the  "finite  part"  integration  techniques 
of  Appendix  B  will  be  applied  when  there  is  no  singularity  present.  On  the  other  hand,  if  t  is 
too  small  then  the  proximity  effects  of  the  approaching  singularity  will  cause  the  integration 
techniques  of  Appendix  A  to  lose  numerical  accuracy. 

Giesing  el  al.  [4]  have  investigated  this  situation.  They  found  that  numerical  difficulties 
arose  in  their  integration  algorithm  when  rmi„/A.r  =  0 - (X)  125.  With  the  present  method,  no 
numerical  difficulties  were  encountered  for  cases  with  rmin/A.v  =  0-01.  Hence,  the  value  of  < 
was  set  at  e  =  0  01,  this  being  more  conservative  than  the  equivalent  value  used  by  Giesing 
el  al.  [4].  Since  A s  is  small,  rmin  =  Are  is  very  small  and  can  be  assumed  to  be  zero;  i.e.  the 
planar  case.  The  small  value  of  e  ensures  that  the  transition  from  the  case  of  very  nearly  coplanar 
surfaces  to  the  case  of  coplanar  surfaces  is  achieved  without  a  discontinuity  in  the  results. 
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6.  COMPARISON  WITH  OTHER  METHODS 


The  method  described  above  has  been  implemented  at  ARL  into  the  FORTRAN  77 
program  DOULAT.  The  formal  of  the  input  data  to  DOULAT,  and  the  definition  ot  pressures 
and  generalised  forces,  is  the  same  as  that  used  in  [lj. 

In  order  to  test  and  verify  the  operation  of  the  present  method,  a  number  of  lifting  surface 
configurations  have  been  investigated.  The  test  problems  selected  cover  a  wide  range  of  frequency 
parameters,  Mach  numbers  and  variations  in  geometry.  These  are  designed  to  test  both  the 
planar  and  nonplanar  aspects  of  the  computational  procedure,  since  both  have  been  significantly 
modified. 

The  various  configurations  are  described  in  Sections  61  to  6-5.  and  the  comparisons  of 
the  results  with  those  of  other  workers  are  shown  in  fables  1  to  7.  The  comparisons  include 
results  obtained  by  a  variety  of  computational  formulations  of  the  doublet  lattice  method, 
lifting  surface  methods,  and  also  some  experimental  results. 

In  the  Tables  the  results  are  presented  in  terms  of  modulus  and  phase.  Where  the  frequency 
parameter  is  zero  ( k  =  0),  the  steady  state  Kernel  function  was  used. 


6,1  AGARD  Horizontal  Wing  and  Tailplane 

The  AGARD  horizontal  wing  and  tail  combination  has  been  specified  by  AGARD  for  the 
calculation  of  generalised  airforces.  The  configuration  is  shown  in  Figure  5.  Calculations  have 
been  made  for  two  values  of  frequency  parameter,  k  =  0  and  k  -  I  •  5.  at  Mach  number  M  =  0-8. 
The  modes  of  oscillation  are  anti-symmetric  and  are  defined  by 
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where  the  origin  of  the  coordinates  is  at  the  apex  of  the  wing.  The  first  mode  is  torsion  of  the 
wing  about  an  axis  38",,  of  the  local  chord,  coupled  with  roll  of  the  tail.  The  second  is  para¬ 
bolic  bending  of  the  wing  coupled  with  pitch  of  the  tail. 

A  number  of  variations  have  been  investigated  for  /;  =  0  0.  0-01.  0  04.  01.  0-2.  0'?.  0-4. 
0-5  and  0-6.  Tables  I  to  3  compare  the  results  obtained  by  different  workers  [2.  5.  8.9,  10] 
with  those  calculated  using  the  present  method.  The  results  of  Davies  [8]  and  Albano.  Perkinson 
and  Rodden  [9]  are  based  on  lifting  surface  methods,  and  the  remainder  are  based  on  various 
computational  formulations  of  the  doublet  lattice  method. 

Using  the  present  method  generalised  forces  have  been  calculated  for  two  different  panel 
distributions.  Pane!  distribution  (I)  consisted  of  8  evenly  spaced  panels  along  the  semi-span 
and  6  evenly  spaced  panels  along  the  chord  of  the  wing,  and  of  8  panels  along  the  semi-span 
and  4  panels  along  the  chord  of  the  tailplane,  also  evenly  spaced.  This  distribution  was  the  same 
as  that  used  by  Albano  and  Rodden  [2],  Panel  distribution  (2)  had  8  evenly  spaced  panels  along 
the  chord  of  the  wing  and  tailplane.  All  panels  on  the  wing  and  tailplane  are  aligned  in  stream- 
wise  strips. 

Table  1  presents  generalised  forces  obtained  for  two  values  ot  frequency  parameter,  k  ----  0 
and  k  =  lx  for  the  coplanar  ease.  The  results  of  the  present  method  agree  well  with  those 
ot  other  workers.  Note  that  the  results  of  Farrell's  doublet  lattice  calculations  were  obtained 
using  an  array  of  It)  panels  along  the  semi-span  and  5  panels  along  the  chord  of  the  wing  and 
the  tailplane. 

It  is  worth  noting  that  in  order  to  converge  to  Hedman’s  [II]  vortex  lattice  results  for  steady 
Mow.  and  to  improve  the  approximation  of  liquation  2-4.  Albano  and  Rodden  found  it  necessary 
to  subtract  the  steady  part  (k  ---  0)  from  the  Kernel  function  A  before  applying  their  integration 
tormula.  and  then  to  add  the  effect  of  a  horseshoe  vortex  which  was  calculated  analytically. 
When  comparing  the  results  of  the  present  method  <  1 )  with  those  of  Albano  and  Rodden  (the 
panel  distributions  for  these  two  cases  being  identical)  for  the  steady  planar  ease  (A  -  0.  h  --  0). 
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it  is  evident  that  there  is  excellent  agreement.  There  is  also  negligible  difference  between  the 
results  of  panel  distribution  (1)  and  (2). 

When  comparing  the  results  of  the  present  method  for  k  =  1-5  and  h  =  0  with  those  of 
the  refined  doublet  lattice  method  of  Giesing,  Kalman  and  Rodden,  it  is  seen  that  results  obtained 
with  the  finer  panel  distribution  of  the  present  method  (2)  are  in  excellent  agreement.  The  results 
of  the  present  method  (1)  do  not  display  such  good  agreement,  although  they  differ  by  a  maxi¬ 
mum  of  less  than  1%. 

It  is  interesting  to  note  that  for  the  unsteady  case  there  is  always  excellent  agreement  in 
the  phase  angle  calculated  by  the  four  formulations  of  the  doublet  lattice  method  for  any  given 
generalised  force.  The  small  variation  that  is  evident  is  very  much  less  than  that  occurring  in 
the  magnitudes  of  the  generalised  forces.  It  is  also  evident  that  the  variation  in  results  obtained 
using  the  present  method  and  panel  distributions  (1)  and  (2)  is  much  greater  for  the  unsteady 
case  than  it  is  for  the  steady  case. 

Table  2  presents  results  for  generalised  airforces  calculated  with  h  =  0-6  for  A:  =  0  and 
k  =  1-5.  The  results  of  the  present  method  once  again  agree  well  with  those  of  other  workers. 

Table  3  gives  a  comparison  of  Davies  lifting  surface  results  for  k  =  1  -5  and  h  =  0  to 
h  =  0-6  with  results  from  the  present  method  obtained  using  panel  distribution  (2).  There  is  good 
agreement  in  the  trends  predicted  for  the  generalised  airforces  as  the  vertical  separation  is  varied, 
together  with  an  acceptable  agreement  between  individual  pairs  of  generalised  forces.  As  described 
in  Section  4,  early  formulations  of  the  doublet  lattice  method  experienced  difficulties  with  non- 
planar  wing  and  tailplane  combinations  where  the  vertical  separation,  h,  was  small  [3.  4],  In  the 
critical  region  of  interest  when  /t  =  0  01  and  h  —  0  04  there  is  acceptable  agreement  between  the 
results  due  to  Davies  and  those  of  the  present  method.  No  numerical  difficulties  were  encountered 
when  the  above  two  h  values  were  used. 


6.2  Stark's  Swept  and  Tapered  T-lail 

The  swept  and  tapered  T-taii  analysed  by  Stark  (12)  is  shown  in  Figure  6.  For  this  /-tail 
the  trailing  edge  of  the  stabiliser  extends  beyond  the  fin  trailing  edge  at  the  stabiliser  fin  junction. 
The  T-tail  is  assumed  to  be  oscillating  in  three  rigid-body  modes.  These  consist  of  yawing  about 
a  vertical  axis  through  the  centre  of  the  root  chord  of  the  tin  (positive  nose  right),  sidesway 


(positive  left)  and  rolling  about 

the 

tin-stabiliser  i 

intersection  (positive  right 

stabiliser  down). 

These  modes  are  defined  as 

/i(.v.  y.  :) 

=  3< 

a  -  0- 15577) 

on 

t  he 

tin 

6-5 

=  0 

on 

t  he 

stabiliser 

fs  ■  (i 

,/>(  v.  y.  ;) 

=  1  ■ 

0 

oil 

the 

tin 

6-7 

=  O' 

0 

on 

the 

stabiliser 

(s  ■  8 

/3(  .v.  l  ,  r) 

—  __ 

.  - 

on 

the 

fin 

6-9 

---  V 

on 

the 

stabiliser 

6-10 

Figure  7  shows  the  panel  distributions  used  on  the  fin  and  stabiliser.  The  present  idealisation 
is  similar  to  that  used  by  Kalman.  Rodden  and  Giesing  [101. 

Table  4  presents  generalised  airforces  calculated  for  A  -  0  and  M  =  0  and  A I  -  0-8. 
The  results  of  the  present  method  agree  well  with  those  of  other  workers.  When  comparting  the 
results  due  to  Kalman.  Rodden  and  Giesing  [!(>|  with  those  of  the  present  method,  the  largest 
difference  in  the  generalised  forces  occurs  for  the  Q , ,  term.  This  term  represents  the  tin  yawing 
moment  due  to  yaw.  and  since  the  side  force  acting  on  the  fin.  (G,.  agrees  well,  this  indicates 
that  the  predicted  chordwise  location  of  the  centre  of  pressure  will  differ.  A  similar  effect  occurs 
for  the  (?,,  term  in  the  results  presented  for  k  -  0-6  and  k  0  9  in  Table  5. 

Table  5  presents  generalised  airforces  calculated  for  M  ^0-8  and  k  --  0-6  and  k  -  0-9. 
The  results  due  to  (i)  Davies,  (ii)  Zwaan.  and  (lii)  Kalman.  Rodden  and  Giesing.  were  obtained 
from  [10], 
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The  results  shown  in  Table  5  compare  favourably  with  those  predicted  by  the  other  methods; 
however,  the  agreement  between  the  magnitudes  of  the  generalised  forces  Q, ,,  Ql2  and  Qn 
as  calculated  by  Farrell  is  comparatively  poor. 

The  generalised  forces  Qxl,  QX1  and  Q 13  represent  the  yawing  moment  due  to  yaw,  sides- 
way  and  roll,  respectively.  Nine  and  eleven  panels  were  used  down  the  chord  of  the  fin  and 
stabiliser  in  the  present  method,  compared  with  only  four  panels  down  the  chord  of  each  surface 
for  the  results  due  to  Farrell,  whilst  the  spanwise  distributions  were  similar.  Since  the  present 
method  has  a  significantly  finer  chordwise  panel  distribution,  it  is  to  be  expected  that  the  yawing 
moments  will  be  predicted  with  greater  accuracy.  As  shown  in  [1],  if  the  chordwise  panel  distri¬ 
bution  is  chosen  to  approximate  more  closely  that  used  by  Farrell,  then  the  results  of  estimates 
for  Q ,,,  Qn  and  £)u  approach  those  obtained  by  Farrell.  Note  that  this  trend  is  also  followed 
by  the  remaining  six  generalised  forces. 


6.3  NLR  Clean  Wing 

This  wing  (Fig.  8)  represents  that  of  the  F- 5  aircraft  and  it  has  been  the  subject  of  both 
theoretical  analysis  and  wind  tunnel  testing  at  the  National  Aerospace  Laboratory  (NLR)  in 
the  Netherlands  [Ref.  13,  Part  11].  During  the  wind  tunnel  tests  the  model  was  oscillated  in 
pitch  about  an  axis  parallel  to  the  >'-axis  and  passing  through  the  mid-point  of  the  root  chord. 

Two  combinations  of  Mach  number  and  frequency  were  investigated.  These  are  M  =0-6 
and  F  =  20  Hz  and  also  M  =  0-6  and  F  =  40  Hz.  The  nondimensional  in-wind  vibration 
modes  are  respectively  [Ref.  13.  Part  11] 

/,(.v,  ;■)  =  (0-330  -  1  ■  00 1  v  +  0-002;-  -  O  OI3.v;-  +  0-052;  ’  -  0-067.vr)  /  6-11 

/2(.v,  ;•)  =  (0-312  —  0-954.V  -  0 - 1 48;  t-  0-532.vc  -  0-339;--  +  0-602.V y2)  /  6  12 

where  the  reference  length.  /.  is  taken  to  be  equal  to  the  root  chord  (/  =  0-639  ni).  Note  that  in 
Equations  6-11  and  6- 12  the  v-  and  ;--ordinates  must  be  given  in  metres.  The  expressions  for 
/,  and /2  assume  that  there  is  no  elastic  deformation  in  the  chordwise  direction.  A  parabolic  defor¬ 
mation  is  assumed  to  occur  in  the  spanwise  direction.  The  major  component  of  both  modes 
is  due  to  rigid  body  pitch  (rotation)  about  an  axis  through  the  centre  of  the  root  chord,  with 
only  a  small  proportion  being  due  to  a  spanwise  bending  and  torsion  contribution. 

The  in-wind  vibration  modes  for  the  clean  wing  are  illustrated  in  Figure  9.  The  normalis¬ 
ation  of  the  modes  is  carried  out  such  that  at  the  chordwise  wing  section  containing  the  point 
(a.  ;  )  =  (0-5648.  — 0  0977)  the  chordwise  slopes  of  the  oscillatory  modes  are  given  by 
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Note  that  the  displacements/,  and/,  are  positive  in  the  positive  r-direction. 

The  panel  distribution  used  in  the  present  method  (Fig.  10)  is  the  same  as  that  used  by 
NLR  [Ref.  13,  Part  II].  The  present  method  was  used  to  calculate  the  unsteady  spanwise  normal 
force  and  pitching  moment  distributions.  The  sectional  normal  force  and  pitching  moment 
are  given  by 


*i  =  2  {'<  C'’°c' 
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where  Zt  is  positive  in  the  positive  z-direction  and  is  the  pitching  moment  about  the  quarter- 
chord  point  of  the  section  (positive  nose  down).  The  chord  c,  is  the  local  chord  of  the  section. 
The  sectional  normal  force  and  pitching  moment  coefficients  are  given  by 

„  2/ 2  K 

Cz j  =  - -r-  >  Ajj  Ati  617 

n  cj  dyi  1  J 

i 

CMJ  =  2,  ;'ij  A,i  (  Vii  “  *i)!l  618 

I 

where  Ay  is  the  nondimensional  pressure  acting  on  a  panel  in  the  /“th  row  of  the  /th  section; 

l2Atj  is  the  area  of  the  panel  in  the  /’th  row  of  the  /th  section; 

*ij  is  the  dimensional  A-ordinate  of  the  mid-span  quarter-chord  point  of  the  panel 
in  the  /’th  row  of  the  /th  section; 

•?j  is  the  dimensional  .v-ordinate  of  the  quarter-chord  point  of  the  /th  section; 

dyt  is  the  dimensional  width  of  the  panels  that  comprise  the  /th  section; 
x  ' 

and  the  summation  /  is  carried  out  over  the  /  panels  down  the  chord  of  the  /th  section. 


Figures  11  to  14  show  the  unsteady  normal  force  and  pitching  moment  distributions  for 
the  clean  wing  obtained  by  the  present  method  and  compare  them  with  results  of  NLR  doublet 
lattice  calculations  and  experimental  measurements.  The  results  of  the  present  method  are  in 
excellent  agreement  with  the  results  of  NLR.  Both  methods  predict  forces  that  are  greater  than 
those  obtained  from  NLR  experiments. 


6.4  NLR  Wing  with  Tipstore 

This  configuration  was  studied  theoretically  and  experimentally  at  NLR  [Ref.  13.  Part  111). 
It  comprises  a  tip  store  added  to  the  clean  wing  described  in  the  previous  section.  The  store 
represents  an  AIM-9J  missile  and  launcher.  The  panel  distribution  used  (Fig.  15)  was  chosen 
to  conform  as  closely  as  possible  to  the  panel  distribution  presented  in  (Ref.  13.  Part  III]. 

The  nondimensional  in-wind  vibration  modes  are  as  follows  (Rel.  13.  Part  111] 

/ 1  (  v,  y )  =  (0  336  -  0-997.V  +  0  045,-  +  0  032 at 

+  0 ■  1 4 1 ’  +  0-001  vy2)  /  (on  wing)  6  19 

=  (0-389  —  1  022.\ )  /  (on  store) 

/(  v,  y)  =  (0  -351  -  I  022*  +  0- 1 62. i  -  0  -296.v.r 

+  0-53Iy2  —  0-0677i-2)  /  (on  wing)  6-20 

=  (0-459  -  I  -075a)  /  (on  store) 

where  the  reference  length  /  is  taken  to  be  equal  to  the  root  chord  (/  =  0-639  m).  Mode/, 
corresponds  to  M  =  0-6  and  F  =  20  H/.  and  mode  /  corresponds  to  M  =  0-8  and  V  =  20  Hz. 
Note  that  the  ,v-  and  y-ordinates  are  to  be  given  in  metres  an  I  the  term  “on  store"  (used  in 
Equations  6-19  and  6  - 20)  refers  to  the  complete  store  and  launcher  assembly. 

The  expressions  for  /,  and  /2  assume  a  parabolic  deformation  in  the  spanwise  direction 
and  no  elastic  deformation  in  the  chordwise  direction.  The  major  component  of  both  modes 
is  due  to  rigid  body  pitch  (rotation)  about  an  axis  through  the  centre  of  the  root  chord,  and 
a  small  proportion  is  due  to  spanwise  bending  and  torsion.  The  displacements  /  and  /,  are 
positive  in  the  positive  z-direction. 
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The  modes  are  normalised  such  that  at  the  chordwise  wing  section  containing  the  point 
(x,y)  =  (0  -5648,  -0-0977)  the  slopes  of  the  oscillatory  modes  are  given  by 

-£(//,)  =-l  6  21 

ax 

^(//2)=~  I  6-22 

ax 

The  nodal  lines  for  the  in-wind  vibration  modes  for  the  NLR  wing  with  tipstore  are  pre¬ 
sented  in  Figure  16.  The  figure  shows  that  as  the  Mach  number  is  increased  from  M  =  06 
to  M  =  0-8  the  nodal  line  is  bent  further  rearwards.  Comparing  Figure  16a  with  Figure  9a 
indicates  that  the  wing  with  tipstore  possesses  a  nodal  line  that  is  bent  further  rearwards,  especially 
over  the  outer  part  of  the  wing. 

The  present  method  has  been  used  to  calculate  the  unsteady  spanwise  normal  force  distri¬ 
bution  on  the  wing  in  the  presence  of  aerodynamic  interference  from  the  tipstore.  The  normal 
force  is  defined  in  Equation  6-17.  Figure  17  compares  the  unsteady  local  normal  force  distri¬ 
bution  calculated  by  the  present  method  with  that  obtained  by  NLR  doublet  lattice  calculations 
and  experiment  [Ref.  13,  Part  III],  It  is  seen  that  the  present  method  agrees  very  well  with  the 
NLR  doublet  lattice  results.  When  compared  with  theory,  the  experimental  results  yield  smaller 
values  for  the  real  part  of  Czj  (especially  near  the  wing  tip)  and  larger  values  for  the  imaginary- 
part  of  Cz j.  This  corresponds  to  a  reduction  in  the  magnitude  of  the  force  together  with  a  phase 
shift.  A  comparison  of  Figures  II  and  17  indicates  that  the  tipstore  acts  as  an  endplate,  thus 
increasing  the  load  on  the  wing. 

The  normal  force  coefficient  Cz  and  pitching  moment  coefficient  C\,  have  been  calculated 
for  the  tipstore  (together  with  its  launcher)  using  the  present  method.  The  normal  force  and 
pitching  moment  about  the  balance  centre  (positive  nose  up)  are  defined  by  the  following 
equations 


/.  =  -P  l2<\v  C ,  0 
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K  ,  , 
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where  c  is  the  mean  geometric  chord  (c  =  0-4183)  and  s  is  the  semi  span  (.v  =  0 - 6226).  The 
normal  force  and  pitching  moment  coefficients  are  given  by 
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where  7.k  is  the  nondimensional  pressure  on  the  A'lli  panel: 

>'k  is  the  dihedral  angle  of  the  A’th  panel; 

vk  is  the  dimensional  v-ordinate  of  the  mid-span  quart'  r-chord  point  on  the  A'th  panel : 
vBAL  is  the  dimensional  v-ordinate  of  the  balance  centre  (vBAL  =-  0  4X0  m); 

/2<4k  is  the  dimensional  area  of  the  A'th  panel; 
and  the  summation  is  taken  over  all  panels  on  the  tipstore  and  launcher  assembly. 
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Table  6  presents  the  unsteady  normal  force  and  pitching  moment  coefficients  obtained  by 
the  present  method  and  compares  them  with  theoretical  and  experimental  results  obtained  by 
Tijdeman  et  al.  [Ref.  13,  Part  III].  The  coefficients  have  been  calculated  for  a  frequency  of  oscil¬ 
lation  F  =  20  Hz  and  two  Mach  numbers,  M  =  0-6  and  M  =  0-8.  When  compared  with  the 
NLR  doublet  lattice  method,  the  present  method  yields  coefficients  whose  magnitudes  are  11% 
to  31%  greater,  the  largest  variations  occuring  in  the  moment  coefficients.  However,  in  contrast 
to  this,  the  agreement  exhibited  for  the  phase  angles  is  excellent,  the  differences  being 
less  than  1  - 1  . 

For  the  spanwise  normal  load  distribution,  excellent  agreement  has  been  obtained  between 
the  present  method  and  the  NLR  doublet  lattice  method.  Even  at  the  wing  tip,  where  the  non- 
planar  effects  of  the  tipstore  are  greatest,  the  agreement  is  still  very  good.  However,  when  the 
unsteady  force  and  moment  coefficients  for  the  tipstore  arc  compared,  it  is  found  that  the  force 
and  moment  coefficients  are  consistently  larger  when  predicted  by  the  present  method.  On 
average,  the  force  coefficient  Ct  is  13%  greater,  and  the  moment  coefficient  %  is  28",,  greater, 
the  differences  in  phase  angles  being  negligible. 

When  the  theoretical  force  and  moment  distributions  over  the  wing  are  compared  with 
experimental  results,  it  is  evident  that  the  theoretical  methods  overestimate  the  magnitudes  of 
these  forces.  This  feature  was  also  apparent  for  results  for  the  dean  wing,  described  in  Section 
6-3.  However,  it  is  not  possible  to  isolate  a  similar  trend  in  the  force  coefficients  calculated 
for  the  tipstore. 

The  reader  is  referred  to  Section  7  for  additional  comments  related  to  the  differences 
between  the  theoretical  calculations  of  the  present  method  and  the  NLR  doublet  lattice  method. 


6.5  NLR  Wing  with  Underwing  Store 

This  configuration  investigated  by  NLR  [Ref.  13.  Part  IV]  involves  the  addition  of  a  pylon, 
launcher  and  missile  to  the  clean  wing  described  in  Section  6  -3.  I  he  store  repress nts  an  AIM-9J 
missile,  and  is  represented  by  additional  thin  lifting  surfaces.  The  panel  distribution  is  shown 
in  Figure  18  and  conforms  as  closely  as  possible  to  the  panel  distribution  presented  in  [Ref.  13. 
Pan  IV], 

The  wing  model  was  oscillated  in  pitch  about  a  50"  „  root  chord  avis  at  a  Mach  number 
M  =  0-6  and  frequency  h  -  20  H/.  The  resulting  vibration  mode  on  the  store  consisted  of 
both  vertical  and  lateral  motion.  The  nondimensional  in-wutd  vibration  mode  is  defined  as 
[Ref.  13.  Part  IV] 
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where  the  term  "on  store''  refers  to  the  store,  pylon  and  launihcr  assembly.  I  he  vertical  com¬ 
ponent  of  the  modal  displacement  is  given  by  ls.  for  which  the  displacements  are  positive  in 
the  positive  e-direction  I  he  lateral  component  ol  the  mode  o  given  by  /,  .  tor  which  the  dis¬ 
placements  are  positive  in  the  positive  i-dtrection  I  he  reference  length  /  0-6396  in  and 
the  v-  and  r-ordinales  are  assumed  to  be  given  in  metres 

I  he  expression  for/v  on  the  wing  assumes  a  parabolic  deformation  in  the  spanwise  direction 
and  no  clastic  deformation  in  the  chordwise  direction  I  he  ma|or  component  of /K  is  due  to 
rigid  body  pitch  (rotation)  about  an  avis  through  the  centre  ol  the  root  chord,  with  a  small 
contribution  being  due  to  a  spanwise  bending  anil  torsion  distribution.  I  he  lateral  component 
of  the  displacement  mode.  /, .  is  composed  of  rigid  bodv  yaw  and  suleswav  of  the  underwing 
store  assembly. 
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The  mode  is  normalised  such  that  at  the  chordwise  wing  section  containing  the  point 
(x,  y)  =  (0-5648,  —0  0977)  the  slope  of  the  oscillatory  mode  is  given  by 


6-29 


on  the  wing. 

The  nodal  line  for  the  in-wind  vibration  mode  for  the  NLR  wtng  with  underwing  store  is 
presented  in  Figure  19.  Comparison  of  Figure  9a  with  Figure  1  9  shows  that  the  nodal  line  of 
Figure  19  is  very  similar  to  that  of  the  clean  wing  at  Mach  number  M  =  0-6  and  frequency 
F  =  20  Hz.  Note  that  the  vertical  displacement  of  the  store  was  taken  equal  to  the  displacement 
of  the  wing  at  location  (x,  r)  =  (0-430.  0-477). 

The  unsteady  spanwise  normal  force  distribution  on  the  wing  in  the  presence  of  the  under¬ 
wing  store  has  been  calculated  by  the  present  method.  The  equation  for  the  normal  force  distri¬ 
bution  has  been  defined  in  Equation  6-17.  Figure  20  shows  the  unsteady  local  normal  force 
distribution  calculated  by  the  present  method  and  compares  it  with  the  results  of  NLR  doublet 
lattice  calculations  and  experimental  data  [Ref.  13.  Part  IV].  As  for  the  two  previous  cases 
involving  the  NLR  wing  combination  (described  in  Sections  6-3  and  6-4).  excellent  agreement 
has  been  obtained  between  the  present  method  and  the  NLR  doublet  lattice  method  for  calcul¬ 
ations  of  the  spanwise  normal  load  distribution  on  the  wing  with  underwing  store.  Even  at  the 
discontinuity  caused  by  the  presence  of  the  pylon,  and  where  the  nonplanar  effects  are  the 
greatest,  the  agreement  is  very  good.  Once  again,  the  experimental  results  yield  a  force  distri¬ 
bution  whose  magnitude  is  less  than  that  predicted  by  the  theoretical  methods. 

When  compared  with  the  loading  distribution  on  the  clean  wing,  the  nonplanar  interference 
due  to  the  pylon  and  store  placed  beneath  the  wing  has  a  pronounced  effect.  The  aerodynamic 
interference  is  responsible  for  an  increase  of  the  loading  inboard  of  the  pylon  and  a  decrease 
on  the  outboard  side.  A  discontinuity  or  jump  in  the  real  part  of  the  normal  load  distribution 
also  occurs  at  the  pylon  station.  However,  no  such  discontinuity  is  present  in  the  imaginary  part. 

The  forces  and  moments  acting  on  the  store  and  pylon  assembly  have  been  calculated  using 
the  present  method,  and  are  defined  by  the  following  equations 
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Cf  is  the  normal  force  coefficient  (positive  in  the  positive  r-direciioni. 

CM  is  the  pitching  moment  coefficient  about  the  balance  centre  (positive  nose  up » : 

is  (he  side  force  coefficient  (positive  in  the  positive  i -direction I ; 

is  the  yawing  moment  coefficient  about  the  balance  centre  (positive  nose  in  the 
positive  e-direction); 

<"  is  the  mean  geometric  chord  (f  -  0  4183  ml 
v  is  the  semi-span  (  v  =  0-h226  m). 
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The  force  and  moment  coefficients  are  given  by: 


„  21 2  \  . 

C/  —  /a  cos  yk  Ak 

res  /  i 

k 


41'  N. 

C m  =  -  —  N  'i  o 


ire's 


cos  yk  Ak  (  vk  -  vBAI  I  / 


2I1  \  . 

Cy  - - 7-  N  M  sin  yk  ,4k 

ItC.V  . 

k 


6*34 


ft -35 


ft  ■  36 


rN 


4/'  \ 

-rr  '-k  sin  Vk  ''k  <  vk Mi\i  l  / 
it  c "  ,v  , 

A 


where  vBAL  is  the  dimensional  v-ordinate  of  the  balance  centre  ( vHA1  ---  0-430  ml.  The  variables 
/.k,  yk,  Ak  and  ,vk  are  as  defined  in  Section  6-4  (liquations  6-25  and  6-26).  and  the  summation 
is  taken  over  all  panels  on  the  store,  launcher  and  pylon  assembly. 

Table  7  presents  the  unsteady  force  and  moment  coetlieients  obtained  by  the  present  method 
and  compares  them  with  theoretical  and  experimental  results  obtained  by  Tijdeman  cl  ill.  [Ref.  1 3. 
Part  IV],  When  compared  with  the  NLR  doublet  lattice  method,  the  present  ntethou  yields 
coefficients  which  are  from  I  -6"„  to  46- 1",,  greater  in  magnitude.  In  contrast  to  this,  the  phase 
angles  exhibit  considerably  better  agreement,  the  differences  being  less  than  2  - 05  .  This  behaviour 
is  similar  to  that  described  for  the  tipstore  in  Section  6-4. 

for  a  further  discussion  of  the  differences  between  theoretical  calculations  with  the  present 
method  and  the  NLR  doublet  lattice  method,  the  reader  is  referred  to  Section  7. 


6.6  Aligned  Panels  for  NI.R  V\ing  with  Store  Configurations 

In  the  application  of  the  doublet  lattice  method  to  lifting  surface  combinations  that  arc 
coplanar  it  o  essential  that  the  panels  be  aligned  in  strc.imwi'c  sirips.  I  or  nonplanar  configur¬ 
ations  this  requirement  may  be  relaxed  as  long  as  the  wakes  ot  panels  do  not  cross  and  cause 
a  singularity  that  is  not  located  at  the  centre  of  the  integration  interval  of  the  panels.  It  must, 
however,  be  toted  that  for  nonplanar  panel  combinations  wh.ch  are  m  close  proximity  to  each 
other,  such  acontigurat ion  will  lead  to  a  more  complex  variation  -  -f  the  planar  and  nonplanar  parts 
of  the  Kernel  function. 

From  an  inspection  of  figure  15  it  is  apparent  that  the  streamw'C  edges  of  panels  on  the 
canard  fins  and  aft  wings  of  the  tipstore  do  not  line  up  n,  .ireumwisk  strips  with  the  panels 
that  comprise  the  launcher  and  missile  body  \  similar  situation  exi'i'  in  I  igare  IS.  where  the 
panels  on  the  canard  tins  and  aft  wings  o!  the  undcrwmg  stoic  do  in  ;  hue  up  m  -tieamwise 
strips  with  the  nearby  panels  on  the  launchei  the  store  bod.  and  the  urns. 

Two  new  panel  distributions  were  created  lor  the  upsi,.  .  amt  pv j. >n  stoic  conltgur.it ions, 
and  they  specifically  embodied  the  featuie  that  all  plana’  and  nonplanar  panels  were  al.gtlcd 
in  streamwise  strips.  I  he  aligned  panel  distribution  loi  a.  .,  g  -.s . r . i  lipstore  is  illustrated  in 
Figure  21.  and  was  created  simplv  In  modifying  the  'Pana'-.  ill  - .  i  ibiit  ion  ol  panels  on  flu- 
launcher  and  missile  bodv,  increasing  the  total  number  o’  p.tiuls  trom  2ss  n>  t(if>.  figure  22 
shows  the  aligned  panel  distribution  lor  the  wing  with  undcrwi  g  store  I  his  new  panel  J ist ri- 
bution  increased  the  number  of  panels  from  244  to  4 i 4 .  and  most  ot  the  retmement  in  the  panel 
distribution  occurred  on  the  launcher,  missile  hodv  and  the  we  e  m  the  vicimlv  of  the  vving- 
pylon  junction. 

Table  K  presents  the  unsteadv  normal  force  and  pud  mg  moment  eoelheients  obtained  for 
the  aligned  panel  distribution  on  the  Nt  R  wing  with  tipstore  using  ilu-  present  method  1  he 
results  arc  compared  with  those  ol  the  M  R  doublet  lattice  method  |f'l  and  earlier  calculations 
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using  the  present  method,  both  of  which  did  not  have  an  aligned  panel  distribution.  The  coeffi¬ 
cients  have  been  calculated  for  a  Mach  number  M  =  0  6  and  frequency  of  oscillation  F  =  20  Hz. 
When  comparing  the  two  results  obtained  by  the  present  method,  it  is  evident  that  the  normal 
force  coefficients  are  in  excellent  agreement.  The  greatest  relative  difference  occurs  between 
the  magnitudes  of  the  pitching  moment  coefficients.  Since  the  normal  force  coefficient  was  less 
sensitive  to  the  change  to  an  aligned  panel  distribution,  the  change  in  pitching  moment  reflects 
a  change  in  the  predicted  centre  of  pressure.  Note  that  it  is  common  to  see  a  greater  change 
in  a  moment  coefficient  than  the  corresponding  force  coefficient  when  a  relatively  coarse  panel 
distribution  is  refined. 

An  inspection  of  Figure  21  shows  that  the  panel  distribution  on  the  surface  representing 
the  launcher  and  missile  body  has  been  considerably  refined.  In  comparison  to  the  non-aligned 
panel  distribution  shown  in  Figure  15.  the  number  of  panels  on  the  launcher  and  missile  body 
has  been  increased  by  a  factor  of  2-5.  Since  this  aerodynamic  surface  is  large  in  relation  lo 
the  complete  tipstore.  and  its  panels  make  a  large  contribution  to  the  pitching  moment,  it  is 
reasonable  to  see  a  large  variation  in  the  pitching  moment  coefficient. 

As  the  aligned  panel  distribution  resulted  in  a  considerably  more  refined  panel  distribution 
for  the  tipstore  assembly,  it  is  not  possible  to  separate  the  effect  of  aligning  the  panels  in  stream- 
wise  strips  from  the  overall  effect  of  having  a  more  refined  panel  distribution.  As  the  results 
obtained  using  the  aligned  panel  distribution  are  similar  to  those  that  would  have  been  obtained 
with  only  a  refinement  in  the  panel  distribution,  there  is  no  evidence  of  any  numerical  difficulties 
associated  with  the  non-aligned  panel  distribution. 

Table  9  presents  the  unsteady  normal  force,  pitching  moment,  side  force  and  yawing 
moment  coefficients  obtained  for  the  NLR  wing  with  underwing  store  using  the  present  method 
and  an  aligned  panel  distribution.  The  results  are  compared  with  those  of  the  NLR  doublet 
lattice  method  [13]  and  earlier  calculations  using  the  present  method,  both  of  which  did  not 
have  an  aligned  panel  distribution.  The  coefficients  have  been  calculated  for  a  Mach  number 
M  -  0-6  and  frequency  of  oscillation  /  =  20  Hz. 

As  for  the  aligned  panel  distribution  for  the  tipstore.  the  aligned  panel  distribution  for  the 
underwing  store  yields  forces  and  moments  that  arc  in  excellent  agreement  with  those  obtained 
for  the  non-aligned  case.  The  largest  relative  differences  occur  in  the  results  for  the  yawing 
moment  coefficient,  and  all  the  phase  angles  experience  only  minor  variation. 

As  mentioned  previously  for  t he  NLR  tipstore.  the  aligned  panel  distribution  is  much 
more  refined  than  the  original  non-aligned  distribution  (see  Fig.  21).  In  particular,  the  number 
of  panels  on  the  launcher  has  increased  by  a  factor  of  three.  Since  this  aerodynamic  surface 
has  a  large  surface  area,  it  is  not  surprising  to  see  that  the  predicted  yawing  moment  coefficients 
have  varied  from  those  obi.,  ned  for  the  non-aligned  panel  distribution.  I  he  relatively  much 
smaller  variation  in  the  pitching  moment  coefficient  is  due  to  the  reduced  level  of  refinement 
of  the  panels  on  the  underwing  store  that  contribute  to  the  pitching  moment  coefficient.  As 
for  the  tipstore  case,  there  is  no  evidence  of  any  numerical  instability  u  ihe  calculations  per¬ 
formed  with  the  non-aligned  panel  distribution 


7.  DISCI  SSION 

Ihe  present  computational  formulation  ol  the  doublet  lattice  method  bus  been  applied  to 
many  dilferent  lilting  surface  configurations.  Ihe  various  combinations  ..nuh-cd  comprised  a 
wing  and  tailplanc  m  tandem,  a  /-tail,  a  wing  in  isolation,  a  wing  with  a  upsiore  and  a  wing 
with  an  underwing  store  ( )t  these  main  configurations  the  latter  two  represent  ihe  most  complex 
notiplanar  aerodynamic  interaction  analysed  here,  and  llicv  'orm  pal  in.u.arly  Jitheiill  lest 
eases  used  m  the  evaluation  of  the  present  method. 

Ihe  wide  range  of  configurations  that  was  analysed  shows  [ h.it  t!u  computer  plogram 
I)()f  I  A  I  can  handle  most  configurations,  provided  that  llicv  are  idealised  with  n  the  limi¬ 
tations  imposed  hy  ihe  d-  uhlcl  lattice  method  Although  not  considered  in  this  report,  full-  or 
partial-span  control  surfaces  may  he  included.  Also,  problems  involving  I  tails  or  twin  tins 
may  be  analysed  without  additional  difficulty. 


When  compared  with  [1J,  the  present  method  has  resulted  in  a  significant  improvement  in 
the  speed  of  computation  of  the  matrix  of  aerodynamic  influence  coefficients.  For  the  NLR 
wing  with  underwing  store  configuration  (described  in  Section  6-5),  the  present  method  reduced 
the  time  spent  in  this  phase  of  the  computations  by  39“-0.  This  improvement  in  computational 
speed  has  been  achieved  without  any  reduction  in  numerical  accuracy. 

The  results  obtained  by  the  present  method  for  planar,  nonplanar  and  intersecting  lifting 
surface  combinations  have  been  compared  with  the  results  of  other  workers  in  the  field.  The 
comparisons  show  that  the  results  of  the  present  method  generally  lie  within  the  range  of  results 
obtained  by  the  other  methods.  Problems  with  a  small  (non-/ero)  vertical  separation  between 
streamwise  columns  of  panels  [3]  can  be  handled  up  to  the  point  where  the  separation  is  negli¬ 
gibly  small. 

However,  for  the  configurations  involving  a  wing  and  tipstore  (Section  6  -1)  and  a  wing 
and  underwing  store  (Section  6-5),  the  comparisons  indicate  that  a  number  of  differences  exist 
between  the  theoretical  results  of  the  present  method  and  those  of  the  NLR  doublet  lattice 
method.  When  compared  with  the  results  of  the  NLR  calculations,  the  present  method  consis¬ 
tently  predicts  forces  on  the  tipstore  and  underwmg  store  that  are  larger  in  magnitude.  There 
does  not,  however,  appear  to  be  any  constant  factor  involved.  Although  the  magnitudes  are 
significantly  different,  the  agreement  in  the  predicted  phase  angles  is  excellent  for  all  the  forces 
and  moments. 

The  possibility  that  the  non-aligned  nature  of  the  pane!  distributions  on  the  NLR  tipstore 
and  underwing  store  configurations  was  contributing  to  the  differences  in  force  and  moment 
coefficients  was  investigated.  Aligned  panel  distributions  for  both  the  tipstore  and  underwing 
store  cases  were  created  (Section  6-6).  and  the  results  of  the  calculations  were  compared  with 
those  obtained  using  the  original  panel  distributions.  The  differences  between  the  two  sets 
of  results  appeared  to  be  due  to  the  refinement  of  the  panel  distributions  which  had  occurred 
for  the  cases  with  the  aligned  panels.  There  was  no  indication  of  any  numerical  instabilities. 

Another  interesting  point  arises  when  the  theoretical  results  for  the  force  and  moment 
distributions  over  the  NL.R  wing  are  compared  with  experimental  results.  It  is  evident  from 
Figures  II.  12.  13.  14.  17  and  20  that  there  is  a  trend  for  the  theoretical  methods  to  predict  force 
and  moment  distributions  whose  magnitudes  are  greater  than  the  experimental  results.  This 
trend  is  generally  followed  by  the  theoretical  results  of  the  present  method  when  predicting  forces 
acting  on  the  tipstore  and  underwmg  store.  A  similar  trend  in  the  NLR  doublet  lattice  predic¬ 
tions  is  not  evident  in  the  results  that  have  been  presented 


8.  COM  I  l  SION 

The  present  formulation  of  the  doublet  lattice  method  has  been  applied  to  a  large  variety 
of  lifting  surlaee  combinations.  As  programmed  at  ARL.  the  method  is  generally  applicable 
to  nonplanar  and  nonparallel  combinations  of  interfering  lifting  surfaces.  I  lie  results  of  the 
present  method  have  been  compared  with  those  obtained  by  other  workers,  and  the  comparisons 
show  generally  good  agreement.  However,  for  the  test  eases  involving  Nl  R  wings  with  either 
a  tipstore  or  an  underwmg  store,  the  agreement  between  the  predicted  magnitudes  of  the  forces 
acting  on  the  store  assemblies  was  not  as  good  although  the  phase  angles  were  in  excellent 
agreement 
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APPENDIX  A 


tutegr  at  .on  pmccdiiie  when  r 1  >  0  over  the  integration  interval 

When  r  ■  0  tin.  ln.m  (  .1: i  integrals  of  Equations  412  and  4-13  no  longer  apply.  Hence, 
the  normalwash  ict'u..  ...  ih-.ient  relating  the  normalwash  on  the  r'th  receiving  panel  to  the 
uniform  procure  a-,  >  ,  i  'lie  »  ill  vending  panel  becomes 


A  I 


where 


A  2 


cos  />, 

4.- 


A,  /, 


A -3 


A  -4 


Ci  -s  .1 

4;: 


A  -  /  ,* 


A  •  5 


IK  :  : 


smaller  snbHitcrv.ds.  v.c  n.av  write 


A  •  (< 


■fr 


where  Ha  line  nn.  ct.d  .  ’  .  s.l  >»er  •/  straight  line  segments  between  points  lm  :  j  \  v,,. 

i  \  \ ,  1  a  1 1  *  I  H .  .  • !'  .mis  I,  ,i  ml  H,  represent  'he  two  v  mis  ol  the  interval  /,. 
I  or  com  eme’ici  i  he  •»  ;  o:-  wll  be  dispensed  with  to  give 


.// 


I  . 

k  : 


,  ji 


A  -S 


A  '> 


Let  the  line  between  points  A:  (vA,  yA.  zA)  and  B  :  (vB, >B,  zB)  be  represented  by  the  para¬ 
metric  relationship 

fab (P)  =  X*b(p)  i  +  KAB(p)  j  +  ZAB(p)  k  -l^p<  I  A  10 

where  *AB(p)  =  *  +  xp  v  =  (vA  +  ,xB)/2  x  =  (rB  —  -vA)/2  All 

*ab(p)  =  >'+>'/>  y  =  (yA  +  >'b)/2  y  =  (yB  -  yA)/ 2  A  - 12 

ZAB(p)  =  2  +  ip  z  =  (zA  +  zB);2  z  =  (zb  -  za);2  A13 

Hence  Equations  A-8  and  A -9  may  be  written  as 


C7(1)(x,(p).  (/>))  d[  d 

y,(p)2  +  z,(p)2  dp  P 


Hi i  =  +  I  G(2>(-V|(p),.V|(p).Z| 
_lj  (y,(p)2  +  z,(/?)2 


z,(/?))  dl 
’)2)2  dP' 


V  l (  Z7 1  —  'r  —  ZAB(  p) 

y.lp)  =  y,  -  yAB(/>) 

4|(/>>  =  A  -  2AB(p) 


dl  I  dr.  a  dr. 


dp  I  dp  dp 


=  (v2  +  f  ’  +  z"2)' 


and  (.r,.  y,,  z,|  is  the  location  of  the  control  point  on  the  r'th  receiving  panel  for  the  combination 
being  considered. 

By  substituting  Equations  A  •  16.  A  ■  17.  and  A  •  IX  into  Equations  A  •  14  and  A  •  1 5  we  obtain 

dl 

=  ~r  — -  dp  A -20 

dp  ,  I  D(p  t 


_  di  " 1  r  6  * 2| 

dp  _  ,  J  D(  / 


Dip)  =  ap ‘  f  bp  t  r 


h  ^  -2|f(>,  -  f)  -t  z(zr  -  f)|  A -24 

<•  =  (y,  —  f>2  +  (zr  —  z)2  A  25 

Now.  in  order  to  evaluate  the  integrals  in  Equations  A-20  and  A  •  2 1 .  let  us  represent  (/' '’(/*) 
and  (j,2'(p)  by  second-order  Lagrangian  interpolation  polynomials 


^  («,/>2  r  />,/>  4  r , )  O' 
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where  the  abscissae  p,  are  llie  /does  oi  (lie  I  .cue  mire  polynomial  I  he  coellicicnts  ht  and  i  , 
mas  bs  calculated  Irom 
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sshere  the  abscissae  />,  aiul  cocltioenis  a..  />  and  .  ,  are  presented  m  I  able  1(1.  By  choosing  the 
abscissae  />,  to  correspond  to  the  abscissae  m  the  common  ( taussian  ipia.lraiare  lormula 


A ;>!.//<  ..  N  : 


It',  'id's; 


eciadon  escuracy  will  he  m.d.ttinsed  eomp;ircil  to  ,:n\  oihd  cln'iee 
•i"-  \  >  •  ■ '  \  y  ,;lo  1  dualio'is  \  Jn  a  I'll  \  2l  ?,yc  oht  i,ri 


3  ■  "*  I  ' oui  1 1|.  t.d'u  ; i ;  uvuei .  ;  u  m 

■:;U ei.ii  o a u s|'oin}mi!  to  the  delude  ;,n. 
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('•  u  d  ""  -h  ■  I/.'  -1,;.  I1 


b 2  —  4  ac  =  0 
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p  ,  1  b  1  , 

-DJp 


p1  p  b  b2  —  lac 

-dp=- - ,  In  I)  +  - ; — 

D  1  a  la2  la2 


—  dp 
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The  indefinite  integrals  corresponding  to  the  definite  integrals  in  Equation  A -38  are: 
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Didp  = 


la  I 


b  +  lap 
(4 ac  —  b2)D  (4«r  —  b2)  I  D 


dp  . 


b2  —  4ac  /  0 


P_d  bp  +  2c _ b _ 

D2  P  (4 ac  -  b2)D  (4 ac  -  b2) 


p-  (b2  —  lac)p  +  be 

-&dp  = - 


2c 


a{4ac  -  h2)D  (4 ac  -  h2)  /) 


—  dp  .  b 2  —  4ac  #  0 


—  dp  ,  b2  —  4ac  ±  0 


I  8  a3 

D1  dp  Mb  +  lap)' 
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t  h  •  lap  r’ 


-  s  dp 
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p 2  _  4a 2  p  b2 

D2  P  (/>  *  lap)2  4 a2 


I) 


dp 


b2  -  4 ac  =  0 


b‘  —  4<;c  =  0 


b2  -  4ai  =  0 


Note  that  since  D2  represents  the  variation  of  r 2  =  r,;  +  r, ;  over  the  interval 
D2  ■  0.  Thus  the  determinant  h2  —  4ac  is  le>s  than  or  equal  to  zero  for  all  practical  i 
interest  and.  therefore,  we  may  neglect  any  integrals  obtained  for  the  case  b2  —  4 ac  - 
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/_  then 
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APPENDIX  B 


Integration  procedure  when  r1  0  in  the  centre  of  the  integration  interval 

In  the  doublet  lattice  method,  when  the  r1  ~  0  singularity  occurs  it  must  do  so  in  the  centre 
of  the  interval  /,  over  which  the  integration  of  Equation  2  •  4  is  defined.  Following  the  method 
of  Appendix  A.  the  interval  /„  may  be  subdivided  into  n  smaller  intervals.  Furthermore,  if  we 
assume  that  the  singularity  falls  symmetrically  within  the  i"th  interval,  then  the  downwash 
influence  coefficient  may  be  written  as 
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77.  A„ 


III 


B  2 


B-3 


The  line  integrals  are  evaluated  over  n  straight  line  segments  between  points  .4m:  (,vAm.  >Am.  trAm) 
and  fi„,:  ( ,vBm.  yBm,  rBm).  Points  A ,  and  B„  represent  the  two  end  points  of  the  interval  /,. 

The  integrals  in  the  first  and  third  terms  of  Equations  B-2and  B-3  may  be  evaluated  using 
the  method  described  in  Appendix  A.  It  now  remains  to  derive  a  method  for  solving  the  finite 
part  integrals  in  Equations  B  2  and  B-3.  If  we  dispense  with  the  subscript  i  then  the  integrals 
of  interest  are 


B  i  (  .,  i  , 

/"'  -  •  - dl 

aJ  r 


114 


dl 


Bo 


Now  since  r  —  0  m  the  centre  of  the  interval  from  point  I :  (  vv  tv  to  point  B  t  i„. 
rH.  rB).  the  line  between  these  two  points  can  be  expressed  by  the  parametric  relationship 
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■  >XI.(/’I 
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--0(2  r 

B  '» 

and  the  control  point  for  the  panel  combination  is  located  tit  (  v  i,.  c,  i.  Vac  that  when  r  1 
the  control  point  on  the  receiving  panel  and  the  lift  point  on  the  sending  panel  (shown  m  I  u:  4 
lie  on  the  same  streamwise  line. 


Using  the  above  definitions  of  XAB(p).  YAB(p)  and  ZAB(p)  we  obtain 


x,{p )  =  -v,  -  (  v  +  .ip) 

y\(p)  =  —>'P 


2|  (/>)  =  -zp 


Hence  we  see  that  r  =  y,2  +  z,2  is  simply 


r2  =  (f2  +  2  V 


and,  as  in  Appendix  A,  we  have  that 


=  (i2  +  p2  +  5 2 ) *  2 


Thus  Equations  B-4  and  B-5  become 


*'!  ayysiJP 

y  +  z  dp  _  |j  p2 


•'  +  z  dp  ~ij  r 


Following  the  method  of  Farrell  [5j,  the  numerators  ( i'"(p I  and  Gi2\p)  may  be  approxi¬ 
mated  by  polynomials  of  order  (//  —  1 1 


G<lV)  =  ^  ^  h„P'  '  G' ”(/>,) 


Gl2\p)  -  \  ^  h„p>  1 


where  the  /jJt  are  the  coefficients  of  the  polynomial.  The  abscissae  />,  are  the  roots  of  the  n’th 
order  Tchebycheff  polynomial  of  the  first  kind.  The  abscissae  may  be  determined  from  the 
following  equation  [Ref.  15.  p,  889] 


( 2/  —  I  )7T 

P,  =  ~  e<>s — ^ - 


By  substituting  Equations  B - 1 7  and  B - 1 8  into  liquations  B - 15  and  B  ib  \xe  obtain 


/'"  =  — , — - — ry  t  ^  C>"G  "</>,) 
yl  +  z  dp  _ _ , 


p2]  =  ---j  ~  ^  r,'21  o1 2 ’(/.,» 

f  +  2  dp  _ _ , 


H  22 


B  23 


11  (P  -  /»„> 

u  ■=  I 

- ./=!,  2 . n  B-24 

fl 

v  =  1 
v  *  i 

As  in  Reference  5.  an  eighth  order  polynomial  has  been  chosen.  I  his  corresponds  to 
n  -  9.  and  the  abscissae  p,  and  the  weights  C,"1  and  C*'1  hase  been  calculated  and  are  listed 
in  Table  II. 
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TABLE  1 


Values  of  Generalised  Airforces,  £M,  for  the  AGARD  Wing-Tailplane  Configuration  (for  A  =  0) 

as  Obtained  by  Different  Workers 


h  =  0 

M  =  0-8 

Op, 

Davies 

Albano 

Perkinson 

Rodden 

Albano 

Rodden 

Giesing 

Kalman 

Rodden 

barrell 

Present 

Method 

(1) 

Present 

Method 

(2) 

On 

L  0,. 

0  4403 
359  9 

0-4425 
359  9 

0-4554 
359  9 

0-4401 

359  9 

0  4377 
360  0 

0-4551 

360-0 

0  4557 
360  0 

*- 

11 

o 

On 

ZOu 

0  6202 
180-0 

0  6121 
180-0 

0  6655 
180  0 

0-6557 

180  0 

0-6457 

180  0 

0  6654 

180-0 

0  •  f>f)53 

180-0 

On 

.02, 

0-1046 

180-3 

0-1054 

180-3 

0-1107 

180-3 

0  1044 

1 80  3 

0-1049 

180  0 

0-1 109 

180-0 

0  1083 

180-0 

022 

<1022 

0  1759 
180  2 

0-1954 

180-2 

0-2237 
180  2 

0-2126 
180  2 

- 

0-2184 

180-0 

0-2235 

180-0 

0-2261 

180-0 

On 
.  On 

1 -5865 
314  2 

1  6022 
314-4 

1  5496 
311-2 

1  5688 
310-7 

1 -4212 
312  3 

1 -5187 
311-7 

1-5714 

310-8 

1=1-5 

0,2 

.  On 

0-9180 

265-5 

0  8910 
266  3 

0-9081 
267  2 

' 

0-9495 

265  4 

0-8752 
265  0 

0-9006 

267-1 

0-9484 

265-8 

O2, 

.  O2, 

1 -0043 
291  -7 

1  0099 

291  4 

1 -0550 
287-2 

1 -0511 
288  •  5 

0  9507 

289-9 

1 -0250 
288-2 

1 -0654 
289-0 

1 -2845 
294-7 

1 -2386 
294  3 

1  2144 

293  •  7 

1  2719 
292-9 

1  1448 

294  -  1 

1  - 1906 
294-0 

1 ■ 2747 

293  2 

TABLE  2 


Values  of  Generalised  Airforces,  for  the  AGARD  Wing-Taiiplane  Configuration  as  Obtained 

by  Different  Workers 


Alhano 

Present 

Present 
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=  0  6 

Perkinson 

Method 

Method 

M 

=  08 

Davies 

Rodden 

Darrell 

(1) 
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359-8 
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360  0 
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0 • 6402 

0-6312 

0-6661 
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0-6868 

Qu 
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0-2527 

0-2674 
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180-0 
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0- 1958 

0-2093 

0  2117 

<>:. 

180-3 

180  2 

180-0 

180-0 

180-0 

1  1009 

1  12(H) 

0  9780 

1  0475 

1  (1838 

301  3 

301  -3 

299  9 

298  X 

298  3 

C',.- 

1  - 1 542 

1  1  128 

1  0755 

1  1072 

1  1 584 

251  ~ 

251  9 

249  5 

250  6 

250-0 

k 

1  5 

0  9i  ("2 

0  9122 

0  8259 

0  895" 

0-9241 

<A-, 

2  "8  0 

7  .  7 

276  9 

275-1 

2"5  -8 

D_. 

1  386’ 

1  3297 

1  2082 

1  260" 

1  5543 

o,. 

289  2 

28"  2 

286  9 

286-4 

TABLE  3 


Generalised  Airforces  for  the  AGARD  W  ing-Tailplane  Configuration  for  Various  Values  of 
Vertical  Separation,  h ,  Between  Wing  and  Taiiplane 
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1  1 
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.  Q 

12 
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V: 
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Davies 
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Method 
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Davies 

.. 
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(2) 
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Method 

(2) 
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i 

Present 
Method 
<  2 1 

0 

1 -5865 
314  2 

1-5714 

310-8 

0-9180 

265-5 

0-9484 

265-8 

.. 

1 -0043 
291-7 

1 -0654 
289  0 

1 -2845 

294-7 

1  2747 
293  -2 

001 

15519 
313  3 

1-5514 

310-7 

0-9332 

264-0 

0-9524 

264-9 

0-9973 

290-4 

1 -0577 
288-7 

1  2890 
294  0 

1  -2736 
292  8 

0  04 

1 -4773 

311-4 

1 -4782 
309-3 

0-9640 
261  -2 

0-9789 

261-7 

0-9843 

287-7 

1  0350 
286-6 

1 -2992 
292  6 

1 

1  2783 
291  3 

0  1 

1-3730 

308-7 

1 ■ 3606 
306-1 

1 -0040 
258-0 

1 -0304 
257-5 

0-9662 
284  -  3 

1 -0043 
282-8 

1  3151 

291-3 

1  2947 

289  2 

0-2 

1 -2624 
305-9 

1 -2463 
303-0 

1 -0479 
255-2 

I -0780 

254-1 

0 • 9428 
281-4 

0-9726 
279  ■  5 

1  3376 
290-4 

1 -3140 
287-8 

0-3 

1  - 1948 
304-0 

1 ■ 1788 
301-2 

1  0790 
253-7 

1  - 1073 
252-3 

0-9271 

279-8 

0-9521 

277-8 

1 -3559 
290-0 

1  327 1 
28'  -2 

04 

1-1511 
302  -  8 

1  1350 
299-9 

1 • 1024 

252  ■  X 
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279  0 

(>•91X4 

276  •  <X 

1 -3695 
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Comparison  of  Unsteady  Loads  on  the  NLR  Wing 
with  Tipstore  Obtained  Using  Aligned  and  Non- 
aligned  Panel  Distributions 
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with  Underwing  Store  Obtained  Using  Aligned  and 
Non-aligned  Panel  Distributions 
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